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We carry out a first-principles density functional study of tiie interaction between a monatomic 
Pt wire and a CO molecule, comparing the energy of different adsorption configurations (bridge, 
on top, substitutional, and tilted bridge) and discussing the effects of spin-orbit (SO) coupling 
on the electronic structure and on the ballistic conductance of two of these systems (bridge and 
substitutional). We find that, when the wire is unstrained, the bridge configuration is energetically 
favored, while the substitutional geometry becomes possible only after the breaking of the Pt-Pt 
bond next to CO. The interaction can be described by a donation/back-donation process similar to 
that occurring when CO adsorbs on transition-metal surfaces, a picture which remains valid also in 
presence of SO coupling. The ballistic conductance of the (tipless) nanowire is not much reduced 
by the adsorption of the molecule on the bridge and on-top sites, but shows a significant drop in 
the substitutional case. The differences in the electronic structure due to the SO coupling influence 
the transmission only at energies far away from the Fermi level so that fully- and scalar-relativistic 
conductances do not differ significantly. 
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I. INTRODUCTION 

Metallic nanocontacts and nanowires are routinely fab- 
ricated by scanning tunneling microscopy (STM) or me- 
chanically controllable break-junctions {MCBJ)} The 
size of these contacts can be as small as a single atom, 
so that electron transport through them is ballistic and 
quantized. The current is carried by quantum channels 
where electrons are partly transmitted and partly re- 
flected. If perfectly transmitting, each (spin degenerate) 
channel contributes to the conductance with one quan- 
tum unit Go = 2e^/h. A massive contact of section area 
A encompasses a channel number N en A and a conduc- 
tance proportional to NGqi^ 

Experimentally one observes that when the contact 
ends are pulled apart the conductance decreases step- 
wise, showing plateaus followed by sudden jumps. In 
general, the values of the conductance plateaus depend 
on the detailed and uncontrolled contact geometry and 
electronic configuration, but averaging over many con- 
figurations makes it possible to plot conductance his- 
tograms which display peaks at most frequent plateau 
values characteristic of the nanocontact of that partic- 
ular metal. ^ The lowest conductance peak is generally 
of order Go — 2e'^/h, thus attributable to a monatomic 
contact in the likeliest geometry preceding breaking. 

Some metals such as gold, platinum and iridium dis- 
play a tendency to form ultimately thin nanocontacts 
consisting of short tip-suspended monatomic chains, 
which have also been visualized by transmission electron 
microscopj«^i^ and/or identified by the presence of very 
long plateaus in the conductance traces versus pulling 



distance. ^'^ In these metals the lowest conductance his- 
tograms peak(s) thus may generally correspond either to 
a chain or to a single atom contact. In Au the peak 
is found at about 1 Gor'^ corresponding to a single 
(spin degenerate) 6s channel. In Pt, where in addition 
to 6s electrons, 5d electrons contribute to the conduc- 
tance, the peak is usually found between 1.5 Go and 
2.0 Go-^^^^^^ 

It has long been recognized that the adsorption of small 
gas molecules on the nanowire or very near the nanocon- 
tact may change substantially the conductance of the 
system, and thus the position of the peaks in the con- 
ductance histograms. For instance, in Pt nanocontacts 
the presence of H2 causes the suppression of the peak 
at 1.5 Go and the formation of new peaks at 1 Go and 
0.1 Go J^i^'^ In presence of CO instead, the lowest con- 
ductance peak is replaced by two new peaks, a major one 
just above 1 Go and a smaller one at 0.5 Go^^"^ 

Such adsorption-induced modifications are usually ra- 
tionalized by guessing a few possible adsorption geome- 
tries and calculating, via density functional theory, their 
Landuaer-Biittiker ballistic conductancCfii In some cases 
the structural assignment is further supported by com- 
parison of the calculated and the measured frequen- 
cies of the molecular vibrational modes. For instance, 
in the case of II2 the peak at 1 Go could be at- 
tributed to a molecular configuration bridging linearly 
the two Pt tips,^^ while in the case of CO it was shown 
theoretically^^ that a conductance value of 0.5 Go may 
be caused by a "tilted bridge" configuration of the CO 
molecule between two Pt contact atoms. In this case 
however, the peak at 1 Go is still unexplained; besides. 
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several other questions remain open, calling for further 
theoretical investigations. To begin with, it remains to 
be clarified what is the electronic mechanism that binds 
a CO molecule to a monatomic wire. Experimentally, 
it has been suggested that the Bly holder model, often 
used to explain the adsorption of CO on transition metal 
surfaceai^ii^ in analogy with metal-carbonyl systems, 
could explain the relative strength of CO absorption on 
Cu, Ni, Pt and Au nanocontacts^ but ab-initio calcu- 
lations proving this point are still lacking. Moreover, in 
the case of heavy atomic species such as Pt and Au it is 
known that spin-orbit (SO) effects will change the elec- 
tronic structure of nanowires^Si^ (for instance altering 
the nature and number of conducting channels, or lead- 
ing a Pt nanowire to become magnetic at distances much 
shorter than those predicted without the explicit inclu- 
sion of SO coupling^°), but it is still unclear how these 
effects would reflect on the chemical binding of a molecule 
such as CO to the nanocontact, or on the ballistic con- 
ductance of the system. 

In order to address some of these points, we will study, 
via ab-initio density functional theory (DFT), a simpli- 
fied model: a monatomic chain of Pt atoms with one 
CO molecule adsorbed on it. We consider a few possible 
adsorption geometries and discuss their electronic struc- 
ture and ballistic conductance (in the absence of tips). 
The role of SO coupling is investigated by comparing 
the results obtained by scalar-relativistic {SR) and fuUy- 
relativistic {FR) pseudo potentials (PP),^^ the latter in- 
cluding the SO coupling. Magnetic effects, which also 
could play a role in this system^S^ are not addressed here 
and will be the subject of future work, along with the 
additional effects caused by the tips. 

We find that when the wire is unstrained and the 
nanowire atoms are kept fixed at their theoretical equi- 
librium distance, adsorption of CO in the bridge site 
(see Fig. [T] in Sec. |lTl is energetically favored with re- 
spect to absorption on the on-top site. The energy of a 
configuration where CO is joining two Pt atoms (substi- 
tutional configuration) is much higher since the energy 
cost of breaking the Pt-Pt bond is large. Hence, this 
configuration could possibly appear only when the Pt-Pt 
bond close to CO is highly strained and almost broken. 
At intermediate strains a "tilted bridge" configuration 
as proposed in Refs. IT5II24I could moreover be energeti- 
cally favorable. In this work we will focus on the up- 
right bridge and substitutional configurations only, two 
examples with completely different symmetry properties. 
Only occasional data on the structure and energetics of 
the on-top and "tilted bridge" geometries will be reported 
for comparison. More results on the on-top configuration 
are given m Refs.imi 

Our results confirm that the Blyholder model~ ex- 
plains quite well the electronic structure in the bridge 
and on-top geometries. In this model both the HOMO 
{5a) and LUMO {27r*) of the molecule form bonding- 
antibonding pairs with the Pt states, giving rise to the 
so-called "donation/back-donation" process (part of the 
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FIG. 1: Schematic representation of the bridge (a) and sub- 
stitutional (b) geometries with the indication of the distances 
that are allowed to vary during the geometrical optimization. 



5(7 density of states, DOS, moves above the Fermi en- 
ergy, Ep, while part of the 27r* DOS moves below Ep). 
In presence of SO coupling this picture retains its valid- 
ity. The states involved in the donation/back-donation 
process are still identifiable and are in fact found to be en- 
ergetically nearly at the same position as in the SR case. 
Similar conclusions are reached in the substitutional case 
although some features are different from the bridge and 
top cases. 

The ballistic conductance of our (tipless) system, 
where the CO adsorbate acts as the scatterer, is found 
to depend on geometry but is not substantially modi- 
fied by the introduction of SO coupling. The theoretical 
SR (FR) conductance is 3.3 Gq (3.2 Go) in the bridge 
configuration, only slightly reduced in comparison to the 
perfect Pt nanowire (whose theoretical value is 4 Go, see 
Ref. i2l[ ). In the substitutional geometry we find instead 
a conductance of 1.1 Go (0.9 Go), thus much smaller 
than the value of the pristine nanowire. Although tip- 
less results are clearly not directly comparable with ex- 
periments, they are consistent with the calculations of 
Strange et al.— , which show a large reduction of conduc- 
tance when the CO passes from the upright bridge to the 
tilted bridge configuration. 

This paper is organized as follows: in Sec.|ll]we briefly 
describe the methods and some technical details. In 
Sec, mil we discuss the geometry and energetics of the sys- 
tem while in Sec. llVl we describe the electronic structure 
of the bridge and substitutional configurations (Sec. lIV Al 
and Sec. IIVBI respectively), reporting both the SR and 
FR results. The ballistic conductance and the energy- 
dependent transmission are presented in Sec. |Vl Finally, 
in Sec. IVII we will draw some conclusions. 



II. METHOD 

We performed our calculations in the standard frame- 
work of density functional theory^! (DFT) using the 
QUANTUM-ESPRESSO package;?^ We used the local den- 
sity approximation^^ (LDA), with the Perdew and 
Zunger- - parametrization of exchange and correlation 
energy. In addition, optimized geometries and adsorption 
energies have been also calculated within the generalized 
gradient approximation (GGA), using th e functional pro- 
posed bv lPerdew. Burke, and Ernzerho^ (PBE). 



3 



Nuc lei and core electrons have been described by 
IVanderbilt ultrasoft pseudopotentiala^^ ( US-PP) , in their 
scalar-relativistic (SR) and fully-relativistic (FR) forms. 
The latter, which include SO coupling effects, are used 
within a two component spinor wavefunctions scheme as 
described in Ref. Hll. LDA PP's of Pt are reported in 
Ref. m, while both for C and for O we generated !JS- 
PP 45 Within GGA we generated SR US-PP for all the 
atoms needed for this work4^ The Kohn-Sham {KS) or- 
bitals are expanded in a plane wave basis set with a ki- 
netic energy cut-off of 29 Ry (32 Ry) in the SR- (FR-) 
LDA case, while a cut-off of 300 Ry has been used for the 
charge density. In GGA calculations (which have been 
performed only in the SR case) we used cut-offs of 32 Ry 
and 320 Ry. The orbital occupations are broadened us- 
ing the smearing technique of Methfessel and Paxton— 
with a smearing parameter a = 0.01 Ry. 

The infinite nanowire has been simulated in a tetrag- 
onal cell, with the wire axis along the z direction. The 
size of the cell along x and y is d± = 18 a.u. ~ 9.5 A, 
so that the interaction energy between the periodic repli- 
cas of the wire is far below 1 mRy. The chemisorption 
energies -Echcm (see Sec. IIIip have been calculated with 
the same d± and the estimated numerical error is less 
than a tenth of eV. The nanowire with one adsorbed 
CO molecule has been simulated by A^pt platinum atoms 
inside a tetragonal supercell, as shown in Fig. [1] for the 
bridge and substitutional geometries. A similar super- 
cell has been used also for the on-top and tilted bridge 
geometries. We checked that the chemisorption energies 
and optimized distances (reported in Sec. lIII[) are at con- 
vergence with A^pt = 15, while for the PDOS (presented 
in Sec. lIVp we used many more Pt atoms in the supercell. 
In order to reduce the effects of the periodic CO repli- 
cas on the bands, we used iVpt = 50 for the FR PDOS, 
and up to Npt = 105 in the SR case. The Brillouin zone 
(BZ) has been sampled with a uniform mesh of fc-points 
along the z direction and the F point in the perpendicu- 
lar directions. In the total energy calculations we used 91 
/c-points for the clean wire (more than enough to ensure 
a converge within 1 mRy for the total energy), while we 
reduced the number of /c-points of the supercell accord- 
ing to A'pt. The PDOS have been calculated employing 
a uniform mesh of (1050/iVpt) fc-points, since the size of 
the BZ reduces linearly with A'pt. 

The ballistic conductance has been calculated within 
LDA, in the Landauer-Biittiker formalism, evaluating the 
total transmission at the Fermi energy. The transmis- 
sion as a function of energy has bee n obtained using the 
method proposed by IChoi and Ihml t^ recently extended 
to the US-PP scheme, both in the SR case^^ and in the 
FR case^^ and implemented in the PWCDND code (included 
in the QUANTUM-espresso package). Since the complex 
band structure calculation (needed to compute the trans- 
mission) is much less sensitive to the periodic replica ef- 
fects (see Sec. |V] for more details), we used 17 Pt atoms 
in the supercell to calculate the ballistic conductance42i 
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FIG. 2: Optimized energy (with respect to the C and O po- 
sitions) for different geometries as a function of the Pt-Pt 
distance, dpt-pt (see text). In the bridge (subs) geometry, 
indicated with solid (dashed) line, the C and O atoms are 
constrained to be on the x [z) axis, while in the tilted ge- 
ometry (dotted line) they can move on the whole xz plane. 
In the SUBS-carbon (-oxygen) the molecule is attached to the 
Pt atom with the C (O) atom. The optimized energy of the 
TOP geometry (only at dpt-pt = dchain) is indicated with an 
horizontal thick line. The zero of the energy is set as the sum 
of the energies of an isolated CO and an isolated Pt chain 
(with the same A^pt and dpt-pt = dchain). This energy has 
been calculated also for some values of dpt-pt ^ dchain (tri- 
angles), and at a large enough dpt-pt it equals -Bbreak- The 
energy gain of placing substitutionally the CO at an already 
broken Pt-Pt bond is then given by -Esubs and is comparable 
to the -Bchcm of CO at the bridge site of an unbroken wire. 



III. GEOMETRY AND ENERGETICS 

We partially optimized the geometry of the system (see 
Fig. [T]) by moving only the C and O atomS(^ while keep- 
ing the Pt atoms fixed on the z axis. In the bridge and 
substitutional configurations, this optimization has been 
repeated for different values of the distance (dpt-pt) be- 
tween the two Pt atoms next to CO while all the other Pt 
atoms are at the theoretical equilibrium distance in the 
isolated chain (dchain — 2.34 A with LDA, dchain = 2.39 A 
with GGA). In the on-top case, only the configuration 
with all Pt atoms at the theoretical equilibrium distance 
has been considered. The GGA total energy of the op- 
timized configurations is shown as a function of dpt-pt 
in Fig. [51 The zero of the energy has been chosen as 
the sum of the energies of the isolated CO and of the 
isolated Pt chain at equilibrium distance (i.e. without 
broken bonds). In the same plot we report also the sum 
of the energies of the isolated CO and of a Pt chain with 
one bond stretched to dpt-pt- 

In Tab. [J we report the chemisorption energies, -Echcm, 
and geometrical parameters calculated for three selected 
configurations: the on-top geometry with all Pt atoms 
equally spaced by dchain, the bridge geometry at dpt-pt = 
dchain, and the substitutional geometry at the optimal 
dpt-pt- For the bridge and substitutional GGA cases, the 
values of i?chcm reported in the table correspond to the 
energy of the bridge and SUBS-carbon curves in Fig. [21 
at the values of dpt-pt pointed by the arrows. 

If the bond is not stretched (dpt_pt — dchain), the 
bridge configuration is favored with respect to the on- 
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TABLE L Relaxed distances (in A) and chemisorption energies (in eV) for the on-top and bridge geometries at dpt-pt ~ dchain, 
and for the substitutional geometry at the optimum dpt_pt- In the substitutional case, the energy cost for breaking the wire 
(calculated as described in Fig. ^ is also shown. 
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top position, the adsorption energy being about 1 cV 
smaller for the on-top configuration within both LDA and 
GGA. The energy of the bridge as a function of dpt-pt 
has a minimum at a distance slightly longer than dchain 
(dpt-pt = 2.56 A in the GGA and dpt-pt = 2.50 A in 
the LDA, not shown here). The substitutional geome- 
try has an energy minimum at a much longer distance 
(dpt-pt = 5.24 A with GGA and dpt_pt = 5.05 A with 
LDA). In this hyper-stretched configuration the substi- 
tutional geometry is favored with respect to the bridge 
configuration. However, here a "tilted bridge" configura- 
tion, where the CO axis lies on the xz plane but is in not 
aligned with the x axis or with the z axis, still has an 
energy slightly lower than the substitutional minimum. 
While this tilted configuration does not reveal an energy 
minimum with respect to dpt_pt, it is anyway preferred 
to the bridge and substitutional configurations at inter- 
mediate distances (about 3.7 A ^ dpt-Pt ^ 5.3 A with 
GGA and 3.8 A «C dpt-pt =^ 5.1 A with LDA). 

Although assessing correct binding energies of CO on 
close-packed Pt surfaces, as the Pt(lll) surface and vic- 
inals to it, can be difficult in DFT with many of the 
commonly used exchange and correlation functionals, po- 
tentially resulting in a wrong predicted site preference, 
in the nanowire case the calculated energy difference be- 
tween bridge and top sites (about 1 eV) is sensibly larger 
than the potential energy corrugation of that problem- 
atic surface case 0.1 — 0.2 eV). Moreover, experimen- 
tally CO has been shown to bind to the bridge site of Pt 
nanowires formed by Pt dimers deposited on a Ge(OOl) 
surface)^ thus partially supporting our prediction for the 
unsupported nanowire. 

In Tab. |T] we also report the energy cost for breaking 
a Pt-Pt bond; it is of the order of 1.9 eV within GGA 
(as can be additionally inferred from Fig. ^ and about 
2.4 eV within LDA. Note that the chemisorption energy 
of substitutional CO on a broken wire (indicated by £'subs 
in the same figure) is comparable to the chemisorption 
energy in the bridge configuration. Moreover, as can 
be seen from Fig. [2 by comparing the two curves SUBS- 
carbon and SUBS-oxygen at a fixed (ipt,_pt, adsorption 
on the C side is much more favorable than adsorption on 
the O side, a standard and chemically obvious finding."'^ 
(filled diamonds and empty diamonds, respectively). 

In all geometries the optimized C-0 bond distance is 
slightly longer than the value dco,oq calculated for the 



isolated molecule, (dco,cq — 1-13 A and dco,oq — 1-14 A 
within the LDA and GGA, respectively). The calculated 
optimal distance between the C atom and the nearest 
neighbour Pt atoms, dpt-c, is about 1.95 A (1.98 A) 
in the bridge geometry within LDA (GGA); it is longer 
than the C-Pt distance found in the on-top geometry, 
dpt-c — 1-82 A (1.85 A). This increase of dpt-c with 
the coordination number of C is in agreement with ex- 
perimental and theoretical values found for CO adsorbed 
on Pt(lll) surfaces. For instance, LEED experiments 
report dpt-c — 1-85 A for atop CO and dpt.^c — 2.08 A 
for CO at the bridge site and DFT calculations at the 
GGA level are quite close to these values (see Ref. [39l 
and references therein). In the substitutional case dpt-c 
is quite similar to the on-top case, while the C-O bond is 
slightly longer than in the other two geometries. Finally 
we note that in LDA bonding distances and chemisorp- 
tion energies do not change significantly in presence of 
SO coupling. This result is in line with the general ob- 
servation, supported by results in the next Section, that 
while Fermi-level related properties including transport 
may be heavily and directly influenced by SO coupling 
through band splittings, the overall energetics, involv- 
ing charge distributions obtained by integration over the 
complex of these bands rather than just those at Ep, is 
much less altered by SO. 



IV. ELECTRONIC STRUCTURE 

In this section we present the electronic structure of the 
bridge and substitutional geometries (calculated within 
LDA), for the specific configurations pointed by the ar- 
rows in Fig. O (dpt-pt = dchain = 2.34 A and dpt_pt = 
5.05 A, respectively). The substitutional geometry pre- 
serves the rotational symmetry of the infinite nanowire, 
and that allows a classification of SR (FR) states accord- 
ing to the quantum number m {rrij) — the projection 
of the orbital (total) angular momentum along the wire 
axis. In the bridge geometry, where the symmetry point 
group reduces to C2v, this is no longer the case. The 
C2v symmetry group contains the elements illustrated in 
Fig. [1] (the identity E, the rotation of tt radians about 
the X axis C2x and the mirror symmetries ai and (72 ), 
but no rotation along z. We analyze the density of states 
{PDOS) projected onto the atomic orbitals of C, of O and 
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FIG. 3: (Color online) PDOS for the SR-LDA bridge config- 
uration. In the upper panel we project onto the even atomic 
orbitals centered on C (solid red lines) and O (dashed blue 
lines) atoms. In the lower panel the projections are onto the 
even atomic orbitals of the Pt atoms below the molecule (solid 
lines). The PDOS for the isolated wire are also shown (dashed 
lines). 



of the Pt atoms in contact with CO, separated accord- 
ing to their symmetry. A similar analysis for the on-top 
geometry is reported in Refs. i25n26i . 



A. CO adsorbed on the bridge site 

We start by discussing the SR electronic structure of 
the Pt nanowire with CO adsorbed at the bridge site 
when all Pt atoms are equally spaced. For arbitrary 
the small group of kz is Cg = {E, cri}, thus we can divide 
the states into "even" and "odd" with respect to the xz 
mirror plane (tri symmetry in Fig. [1]). 

We report in Fig. [3] the electronic PDOS projected on 
the even atomic orbitals of C and O and on the even 
atomic orbitals of a Pt atom in contact with CO. The 
PDOS are decomposed into the projections on the s, px 




FIG. 4: (Color online) Two dimensional contour plot of the 
charge density of bat (left) and 5aa (right) ifS-eigenstates 
calculated at kz — 0. Dots mark the positions of atoms, while 
solid red (dashed blue) lines are positive (negative) isolevels 
with the following values of density: ±0.0005 ±0.001, ±0.002, 
±0.005, ±0.01, ±0.02, ±0.05 electrons/(a.u.)^ 



or Pz orbitals of C and O and into those on the s, d^z^-r^ j 
dj-z or dx2_y2 orbitals of Pt. Some of the peaks corre- 
spond to the even levels of the molecule. The 3a and 
4(7 molecular states are in the C and O s and Px PDOS, 
at —22.5 eV (not shown) and at —10.4 eV, respectively, 
while peaks corresponding to the Itt and to the 27r* states 
are at about —6.3 eV and at 3.3 eV in the C and O Pz 
PDOS. The 4cr state is slightly hybridized with the even 
Pt orbitals, while the 3(t state is not. 

The broadening of molecular features and the presence 
of several additional peaks indicate a chemical interaction 
between the molecule and the wire. Interaction between 
the (7 orbitals and Pt states is visible in the projections 
on the s and Px orbitals of C and O from —7.5 eV up 
to 3 eV. In two energy regions, at about 2.1 eV and 
between —7.5 eV and —6.0 eV there are several large 
peaks in the s projection and in both the s and Px pro- 
jections, respectively, whereas at other energies the pro- 
jections are smaller while not completely vanishing. The 
peak centered at about 2.1 eV in the PDOS projected on 
the C s orbital is due to the 5a orbital and shows that 
there is a depopulation of the molecular HOMO level 
(donation). These states are antibonding Pt-CO states 
while the corresponding bonding states show up as the 
double-peak feature below Ep, between about —7.5 eV 
and —6.0 eV. Thus the 5a orbital is broadened and gen- 
erates a set of occupied bonding states 5ab below Ep and 
a set of empty antibonding states 5aa above Ep. Some of 
the KS-eigenstates associated to these peaks are shown in 
Fig. m In this figure we draw a contour plot of the charge 
density of a Bat state at —6.1 eV and of a 5aa state at 
2.1 eV in the xz plane (which contains both the molecule 
and the wire). We note that the interaction brings part 
of the 5a orbitals below the Itt level, and that there is 
an orbital mixing between Aa and 5a levels. As a conse- 
quence, the projections of 5ab states on the s and Px or- 
bitals of O are higher than those on the s and px orbitals 
of C. These are the same features which characterize the 
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FIG. 5: (Color online) Two dimensional contour plot of the 
charge density of KS-eigenstates (calculated at = 0) cor- 
responding to the bonding (left) and antibonding (right) be- 
tween even 2n* and even Pt states. 
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absorption of CO on transition metal surfaces, well de- 
scribed in the literature (see Ref . ITs ) . The hybridization 
of the 5tT orbital is mainly with the m ~ bands of 
the Pt wire, which have predominant d^^2_r2 character 
at low energies and s character at energies higher than 
0.4 eV. Moreover a smaller, but relevant, hybridization 
of 5(7 states with dxz {\iti\ = 1) Pt orbitals is present both 
in 5i7f, and in Sctq states. 

Broad features and additional peaks due to the hy- 
bridization are present also in the PDOS. Moreover 
some hybridization of tt levels with Pt states occurs also 
at —6.3 eV, which is the position of the Itt states. The 
strong peak at about —2 eV and the weaker features 
above and below —2 eV are due to bonding hybridization 
between the even 2tt* orbital and the Pt states. The pres- 
ence of these new states, which we may call 27t^, shows 
that the 27r'* orbital of CO is partially occupied corre- 
sponding to back-donation. The 27r'^ level broadens be- 
tween 3.2 eV and 3.6 eV and forms the antibonding 27r^ 
states. A contour plot of the charge density associated to 
the 27r'^-derived states is shown in Fig. O a bonding or- 
bital at —2.0 eV on the left and an antibonding orbital at 
3.5 eV on the right. The 27rJ states have a large contribu- 
tion from oxygen orbitals, with very low projections 
on the carbon atom (see also Fig. [3]). This is caused 
by a hybridization between the Itt and 2tt* orbitals of 
the molecule due to its interaction with the wire, and 
has been well characterized for CO adsorbed on metal 
surfaces.™ The 27t^ hybridization at —2 eV involves both 
TO = (s and d^z^-r^) and |to| = 2 (c?3,2_j^2) Pt states, 
and because of this interaction a relevant portion of the 
dx2-y2 PDOS moves outside the energy range of the wire 
|to| = 2 band. The |to| = 1 (d^z) states instead give a 
much smaller contribution at —2 eV and interact in a 
wider energy range between —4 eV and —2 eV. 

In Fig. [5] we show the PDOS projected on the odd 
atomic orbitals. As found in the even tt interaction, in 
the PDOS projected onto the Py orbitals of C and O two 
peaks correspond to the odd Itt molecular level at —6 eV 
and to the odd 27r^ state at 2.8 eV, while additional peaks 



FIG. 6: (Color online) PDOS for the SR-LDA bridge configu- 
ration. The projections are onto the odd atomic orbitals cen- 
tered on C (solid red lines) and O (dashed blue lines) atoms 
and on the Pt atoms below the molecule (solid line, lower 
panel). The PDOS for the isolated wire are shown with a 
dashed line. 



are caused by interaction. Here we find two new peaks 
below E-p, one at —4.3 eV and the other at —2.2 eV. 
They correspond to the odd 27r^ levels and their splitting 
reflects the different hybridization of the 27t* orbital with 
distinct bands of the wire. The peak at —4.3 eV is present 
only in the PDOS projected on dyz, while the peak at 
—2.2 eV is essentially due to an hybridization with the 
|m| — 2 band, since the corresponding dxy peak is much 
stronger than the one in the dyz PDOS. As noticed in the 
even PDOS, this strong interaction brings some charge 
in the 2tt* molecular state and perturbs the PDOS of the 
wire, especially the |m| = 2 component {d^y in the odd 
case) where a lot of states are now outside the energy 
range of the clean wire |m| = 2 band. 

With addition of SO in the FR case this 
donation/back-donation picture does not change. 
Now states fall into the two irreducible representations 
and P'' of the double group CP and it is not possible 
to separate even and odd states. Since in the absence 
of magnetization time reversal symmetry holds, there 
is a (Kramers) degeneracy between the P^ band at kz 
and the P'* band at —kz. Therefore we can analyze the 
CO-Pt nanowire system by focusing on bands belonging 
to one of the two symmetries, for instance the P^ 
bands. In C or in O there are four FR atomic orbitals 
that transform according to P'^, one derived from the s 
{I ~ 0) state, with {j,jy) ~ {i/2, 1/2) and three derived 
from the p [I — I) states, with {j,jy) = (1/2,-1/2), 
(3/2,-1/2), and (3/2,3/2). 

In Fig. [7] we show the PDOS projected onto these four 
P'^ atomic orbitals. The molecular levels can be eas- 
ily identified. At lower energies two sharp peaks are at 
—22.4 eV (not shown) and at —10.5 eV in all the four 
PDOS and correspond to the 3a and 4cr SR states, re- 
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FIG. 7: (color online) PDOS for the FR-LDA bridge config- 
uration (wire along y, see text). The projections are on the 
F"* symmetry FR atomic orbitals of C (solid red lines) and O 
(dashed blue lines). 



spectively. Two other peaks are close together (—6.3 eV 
and —6.0 eV) and are both present in the PDOS pro- 
jected onto the two states with / = 1 and jy = —1/2, 
but only the low energy peak is evident in the PDOS 
projected onto the jy = 3/2 state. The position of these 
two peaks coincides with the even and odd Itt states ob- 
served in the SR PDOS. The odd Itt state (-6.0 eV) has 
very low projection onto the jy = 3/2 state since the for- 
mer is oriented along the xy direction, while the latter 
is made up of m = 1 orbitals, that are oriented in the 
xz plane (the quantization axis is y). The broad feature 
between —7.3 eV and —6.2 eV, which is present in all 
four PDOS, can be distinguished from the neighbouring 
Itt peaks since it has much more weight on the O atom 
rather than on the C atom, similar to the Serf, SR states 
(see PDOS of Fig. [3]). Therefore we can identify this 
feature as the FR analog of the 5(Tb states, which are in 
the same energy range in the PDOS projected on the s 
and Px orbitals. Above the Fermi energy, we can identify 
the empty baa antibonding states which give rise to the 
broad peak at about 2 eV in the jy — 1/2 PDOS. 

In the energy range between —4.3 eV and —1.1 eV we 
find some features which are more evident in the / = 1 
components of the PDOS, and have very low weight on 
the C orbitals. In the SR case the energy range of the 27rJ 
states goes from —4.4 eV (lowest peak in the py PDOS 
of Fig. [5]) to about — 1 eV (tail of the even 27rJ peak in 
the Pz PDOS of Fig. [3]) and they have very low projec- 
tions on the carbon, thus we can identify these features 
in the FR PDOS as the FR analog of those states. The 
corresponding even and odd 27r* antibonding states are 
responsible for the peaks at 3.2 eV and 2.8 eV, respec- 
tively. Although their splitting is too small to be resolved 



O 



: 1 

A , / 


Pt (C side) 

I, ii 


S 

, , V- 


^ A 










\ , 1, 


CO 






Pz 


Pt(Oside) \\ 

-A . . 


S 


'-i 


1 '^--^---------'^^'^^^ 



-12 -10 -8 -6 -4 -2 2 4 
E - Ep [eV] 



FIG. 8: (Color online) PDOS for the SR-LDA substitutional 
configuration. In the central panel the projections are on 
m = atomic orbitals centered on C (solid red lines) and 
O (dashed blue lines) atoms. In the top (bottom) panel the 
projections on m = atomic orbitals of the Pt atom next to 
the C atom (O atom) are indicated by solid lines, while the 
corresponding PDOS of the isolated wire are superimposed 
with dashed lines. 

in a single PDOS with this value of the smearing, only 
the even peak (such as the even Itt) is evident in the 
jy = 3/2 PDOS, while both are present in the Z = 1, 
jy = —1/2 PDOS and both have much smaller weight in 
the I = component. We can therefore conclude that 
although SO changes the symmetry of the orbitals, the 
mechanism of donation and back-donation still describes 
well the bonding between the molecule and the wire. 



B. Substitutional CO 

In the substitutional configuration (see Fig. [1} the ro- 
tational symmetry is preserved, hence SR states with dif- 
ferent TO cannot hybridize. Therefore we can study the 
interaction by focusing on the PDOS projected on m = 
and |to| = 1 orbitals only (Fig.[8]and Fig.[9l respectively), 
disregarding |to| = 2 states that are not present in the 
molecular levels of CO. The Aa and 5a molecular levels 
give rise to two peaks at —11.0 eV and at —8.5 eV in 
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FIG. 9: (Color online) PDOS for the same system as in Fig.[8l 
but here the projections are on |m| = 1 atomic orbitals. 



the C and O s and Pz PDOS. As in the bridge config- 
uration, the 5(Tb level is lower than the Itt derived level 
(see below). We do not find here isolated peaks above 
E-p which correspond to the antibonding 5aa states, but 
an almost flat plateau which extends in the whole energy 
range of the m — bands of Pt. The hybridization be- 
tween the 5cr CO orbital and the Pt states is different on 
the two sides of the molecule, as can be seen comparing 
the PDOS projected onto the m = orbitals centered on 
the two opposite Pt atoms. At —8.5 eV the 5(J CO or- 
bital is more coupled to the Pt on the C side (especially 
via d^z^-r^ states), and this causes a higher depopulation 
of m = states between —7.3 eV and —2 eV on the Pt 
next to the C atom with respect to the Pt on the O side. 

The PDOS projected onto the |m| = 1 {px or py) or- 
bitals of C and O are reported in Fig. [9] and show two 
peaks due to the tt molecular states, one at —6.8 eV 
(Itt) and another at 2.8 eV (27r*). However, in contrast 
with the bridge case, instead of new intense peaks be- 
low the Fermi energy, we find in the substitutional case 
a plateau which spans the whole energy range of the Pt 
wire |m| = 1 band. 

Although the hybridization with Pt states is differ- 
ent from the bridge geometry, donation/back-donation 
is present here too. In fact a portion of the plateau in 
the PDOS on m = states extends above Ep, while the 
plateau on the \m\ — 1 PDOS lies mainly below E-p. In 
order to have an estimate of the donation wc can con- 
sider the integral (from — oo to Ep) of the PDOS on the 
C and O orbitals forming a states in the molecule (s 
plus Px in the bridge geometry, s plus pz in the substitu- 
tional). This integral gives 5.3 both in the substitutional 
and bridge geometries, while in the isolated molecule it 
gives about 5.9. Estimating in the same way the amount 
of back-donation, by integrating the PDOS on the or- 
bitals which form the tt states {pz plus Py in the bridge 
geometry, px plus Py in the substitutional), we find 4.7 



FIG. 10: (Color onUne) PDOS for the FR-LDA substitutional 
configuration. The projections are on the |mj| — 1/2 and 
\mj\ —3/2 FR atomic orbitals of C (solid red lines) and O 
(dashed blue lines). In the case I — 1, \mj\ — 1/2 the PDOS 
projected onto states of different have been added. 



electrons both in the substitutional and in the bridge ge- 
ometry, to be compared to the value of 3.9 obtained with 
the isolated CO. 

In the FR case we can label states according to the 
total angular momentum m-,-, and hybridization occurs 
only among states with the same rrij . In Fig. [10] we report 
the PDOS projected onto the FR atomic orbitals of C and 

separated according to the values of I and rrij of the 
spin-angle function. The peaks in these figures can be 
easily identified and their position compared with that 
of the SR case. In addition to the molecular a levels 
present only in the PDOS on the \mj\ = 1/2 states at 
-23.7 eV (not shown), -11.2 eV and -8.6 eV, there 
are two peaks in the I = 1, \mj\ — 1/2 PDOS, one at 
—7 eV and another at 2.4 eV. They can be matched 
with the two peaks at —6.9 eV and at 2.8 eV in the 

1 — 1, \mj\ = 3/2 PDOS, and correspond to the SO split 
Itt and 2tt* states of the molecule. With respect to the 
isolated molecule, the 1/2-3/2 splitting of the tt states is 
enhanced by the interaction with the Pt states, especially 
for the 2tt* states. 

As in the SR case, the interaction between the molecule 
and the wire is visible here as a plateau which extends 
in a wide energy range. In the PDOS projected onto the 
I — 0, \mj\ ~ 1/2 states the plateau extends between 
—7 eV and 3 eV, and is due to the hybridization be- 
tween the 5cr orbital with the \mj \ = 1/2 bands of the Pt 
nanowire. The small peak at about —3.3 eV and the gap 
just above it are due to an anticrossing of the \mj \ = 1/2 
bands (see Ref. 21), which generates new peaks and gaps 
in the FR-DOS of the Pt nanowire. In the PDOS pro- 
jected on the I = 1, \mj \ = 3/2 states there is a plateau 
between —4.2 eV and 0.8 eV due to the hybridization 
between 7r-derived states and \mj\ = 3/2 states of Pt. 
This plateau is localized in the same energy range as 
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that of the SR TT-\m\ = 1 bands, shghtly increased by the 
spin-orbit spUtting of the |m| — 1 band. A hybridization 
gap of the \mj\ = 3/2 band of Pt is visible just above 
-2 eV. In the PDOS projected on I = 1, \mj\ = 1/2 
states we find contributions from both cr-derived an tt- 
derived states. In summary, also in the substitutional ge- 
ometry we can conclude that the donation/back-donation 
model gives a good description of both the SR and FR 
electronic structure. Near the Fermi level the changes 
caused by SO are not large in extent, but still quite vis- 
ible in the electronic structure. This will reflect on the 
ballistic transport properties, which are discussed in the 
next Section. 



BALLISTIC CONDUCTANCE 




-4 -2 

E - Ep [eV] 



In this section we present the transmission of an in- 
finite ideal platinum chain as a function of energy, for 
the two CO adsorption geometries discussed above, pur- 
posely without tips but with the adsorbed molecule as the 
sole scatterer. This idealized transmission measures the 
amount of obstacle posed by the molecule to electron free 
propagation, in addition to that, molecule-independent, 
caused by the tip-wire contacts - which as was said are 
left out here. 

The transmission is calculated with the method devel- 
oped in Refs. ,21.j35.,40, where the self-consistent poten- 
tials for the left lead, for the scattering region and for 
the right lead are calculated using the geometries shown 
in Fig. [H The self-consistent potential of one unit cell of 
the Pt chain in the leftmost part of the supercell is used 
to calculate the generalized Bloch states of the left and 
right leads, while the rest of the supercell is the scatter- 
ing region. We checked that the generalized Bloch bands 
with real wavevector of the Pt chain calculated with this 
potential perturbed by CO match those calculated with 
the exact potential within 0.05 eV. 

The SR transmission Tsr as a function of energy for the 
bridge geometry is shown in the upper panel of Fig. fTTI In 
the same plot we display also the FR transmission Tfr for 
selected energies and the number of channels available for 
transport {Nsr in the SR case and NpR in the FR case). 
At the Fermi level four (spin-degenerate) SR channels are 
available, three even and one odd. The |m| = 2 bands 
do not cross Ep, hence they do not contribute to the 
conductance. The odd |m| = 1 band is almost perfectly 
transmitted, while the even |to| = 1 band and the two 
m = bands are partially reflected. The resulting con- 
ductance is GsR = 6.6 e'^/h, to be compared with the 
value 8 e'^/h of the clean wire. The FR conductance, 
Gfr. — 6.5 e'^/h, is quite close to the SR value because, 
although the SO splitting of the |m| — 2 bands brings 
one two-fold degenerate channel (with \mj\ = 5/2) close 
to E-p, this additional channel is poorly transmitted (see 
later) and does not modify the conductance. These val- 
ues of conductance are close to those found in the on-top 
geometry: G = 6.1 e^//i in both SR and FR calculations. 



FIG. 11: (Color online) Electron transmission for the bridge 
geometry. In the top panel the SR (FR) transmission Tsr 
(Tfr) is displayed with solid lines (green circles), while the 
number of channels Nsr (Nfr), which equals the transmission 
of an ideal monatomic chain, is shown with dashed lines (short 
dashed green lines) . In the bottom panel we separate the even 
and odd components of the SR transmission, Tcvon and Todd 
(above and below, respectively), as well as the number of even 
and odd channels, A^'even and iVodd- 



The energy dependent transmission shows instead 
more pronounced differences between SR and FR results. 
In the SR case we can separate the contribution of the 
even and of the odd states (central panel and bottom 
panel of Fig.[TTl respectively), while in the FR case states 
belonging to the or to the r4 representations have the 
same transmission. We discuss first the even contribution 
to the transmission, trying to establish a connection with 
the features in the PDOS projected on the correspond- 
ing atomic orbitals in Fig.[3]of Scc. lIV Al Below —4.3 eV 
there is only one m = channel which is almost perfectly 
transmitted above —5.5 eV, while it is partially reflected 
at lower energies and totally reflected at —6.4 eV. This 
happens both in the SR and in the FR cases. We note 
that this energy corresponds to the position of the even 
Itt states, which are rather localized on the CO and on 
the two neighboring Pt atoms and are responsible for the 
large mismatch between the m ~ PDOS of those two 
Pt atoms and the corresponding PDOS of a Pt atom dis- 
tant from CO (compare solid and dashed lines in the s 
and d;}^2_j.2 PDOS of Fig. [3]) At higher energies, where 
in principle a new even channel (with \m\ = 1) becomes 
available, the even contribution to the transmission grows 
significantly only after —3.3 eV, while below this energy 
only the m = channel is transmitted. Actually, in the 
energy region between —4.2 eV and —2.1 eV, the even 
m| = 1 states of the wire are involved in the hybridiza- 
tion with the 27r* molecular orbital and therefore they 
do not transmit (compare the PDOS projected on the 
dxz orbital with that projected on (i3z2_^2, in Fig. [31 in 
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this energy range, the matching with the PDOS of the Pt 
distant from CO is better in the case of d3z^_r2 orbitals 
rather than d^z orbitals). 

Above —2 eV the 27rJ hybridization has more weight 
on the TO = states of Pt than on \m\ = 1 states, hence 
the rise of transmission between — 2 eV and — 1.8 eV has 
to be ascribed to a better transmission of the |to| = 1 
channel, rather than to the additional to = channel 
that is present from —1.8 eV onwards. Actually both 
TO = channels transmit only partially in this region, 
and together contribute just for about one half the value 



of T,, 



Between —1.6 eV and —0.4 eV also the even 



\m\ = 2 channel becomes available, but the transmission 
grows only slightly and remains much smaller then the 
number of available channels. The bad transmission of 
|to| = 2 states can be easily related to the large mis- 
match between the d^2_y2 PDOS of the two Pt atoms 
next to the molecule and that of a Pt atom distant from 
CO. Between —1.7 eV and —1.5 eV there is also a large 
difference between SR and FR transmission since the an- 
ticrossing of the 1/2 and 3/2 bands removes several chan- 
nels from that energy region decreasing the FR transmis- 
sion. Above —0.4 eV the \m\ —2 channel is not avail- 
able anymore and we find a small drop in transmission, 
which then remains almost constant up to Ep. From 
up to 0.8 eV the transmission decreases regularly, since 
first the \m\ — 1 channel (at 0.4 eV) and then one of 
the two m = channels (at 0.8 eV) disappear. Above 
0.8 eV only the to = channel with predominant s char- 
acter is present: its transmission vanishes completely at 
2.1 eV and it shows a dip about that energy. This energy 
value coincides with the position of the Sctq states (see 
the PDOS projected onto the s orbitals in Fig. [3]). This 
dip is present also in the FR transmission although it is 
slightly shifted towards lower energies. 

The odd channels are available and can contribute to 
the transmission between —4.3 eV and 0.8 eV. The 
\m\ = 1 channel has good transmission in the whole 
energy range except for two narrow dips near —2 eV. 
The lowest dip corresponds to the odd 27rJ peak present 
at that energy, while the position of the highest can be 
matched with the energy of a strong peak in the PDOS 
projected on the dxy states of Pt (which does not corre- 
spond to any peak in the PDOS of C or O). As for the 
even |to| = 2 channel, also the odd |to| = 2 channel is 
almost blocked because the corresponding Pt states are 
more perturbed by the interaction with the 2tt* orbitals 
of CO. This difference in the transmission properties be- 
tween |to| = 1 and |to| = 2 odd channels could be pre- 
dicted by the PDOS reported in Fig. [BJ in fact, the dyz 
PDOS looks similar to that of a Pt atom of the pristine 
nanowire, while the dxy PDOS does not match at all the 
original PDOS of the wire. 

A feature similar to the dip at 2.1 eV has already been 
observed in the transmission curve calculated for gold 
nanowires with CO absorbed on-topf^i^ but in that case 
the dip is closer to the Fermi energy. This causes a drop 
in the predicted conductance of the Au nanowire when 
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FIG. 12; (Color online) Electron transmission for the substi- 
tutional geometry. In the top panel the SR (FR) transmission 
TsR (Tfr) is displayed with solid lines (green circles), while 
the number of channels Nsr (Nfr), is shown with dashed lines 
(short dashed green lines). In the middle (bottom) panel we 
separate the SR (FR) transmission and number of channels 
according to the angular momentum \m\ {\mj\). The contri- 
bution from |m| = 2 (|mj| = 5/2) channels to the total SR 
(FR) trasmission is practically zero, hence it is not shown (see 
text). 



CO is adsorbed because, in that case, the 6s state is 
the only conducting channel left at Ep. Instead, in the 
pristine Pt nanowire three more channels are available 
for conduction, and the dip in the to = transmission 
caused by CO adsorption lies above the Fermi level. 

The SR and FR total transmissions as a function of 
energy for the substitutional geometry are shown in the 
top panel of Fig. [T^l In this geometry we find a conduc- 
tance of about 2.2 e^/h, much smaller than in the bridge 
geometry. Actually at the Fermi level the transmission 
of the two degenerate |to| = 1 channels is about one half 
and that of the to = channels is quite small (slightly 
above 0.1). In the other energy regions the transmission 
is even lower, since the to = channels are almost totally 
blocked and the transmission of the \m\ = 1 channels 
remains always below one half. The 5(T-derived states 
could in principle transmit because they hybridize with 
Pt states, but actually this is not the case since the cou- 
pling of the left Pt to CO is quite different from that 
of right Pt, as noted before in Sec. IIVBI Moreover, the 
transmission due to the |to| = 2 bands of Pt (not shown 
here) is close to zero, since the CO has no states with 
matching symmetry. 

Thus the SR and FR calculations give similar results 
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also for the substitutional geometry; the calculated FR 
value of the conductance is 1.7 e'^/h to be compared with 
G ~ 6.5 e'^/h for the upright bridge. In the substitu- 
tional case the contributions from the \mj\ = 1/2 and 
the \mj\ =3/2 channels are almost equal, in agreement 
with the fact that in the SR case only the \m\ — 1 chan- 
nels contribute to transmission. 

We can pinpoint some dips specific to the FR energy 
dependent transmission, that are due to the SO induced 
splittings of the bands. They are evident if we look sepa- 
rately at the \mj\ = 1/2 and \mj \ — 3/2 contributions to 
the total number of available channels: the former goes 
to zero between —3.4 eV and —2.5 eV, while the latter 
vanishes between —2.2 eV and —1.4 eV. The \mj \ = 5/2 
contribution to the transmission (not shown) is almost 
zero, since the \mj\ = 5/2 channels have |to| = 2 orbital 
components that are completely blocked by the CO. It 
should be pointed out here that the details of all channels 
will be modified by the onset of magnetism, ^■^•''■^ and this 
could in principle modify the ballistic conductance of our 
system. This aspect however is beyond the scopes of this 
work. 

Although as was said above our conductance values 
cannot be directly compared with experimental data be- 
cause they do not include the effect of the tips, they do 
indicate that in the substitutional configuration the con- 
ductance is about 3 times lower than in the upright bridge 
configurationfi^ where at difference with the substitu- 
tional the m = channels are not so strongly reflected. 
The Pt states with angular momentum |m| = 2 are (more 
or less strongly) blocked in both cases, but for different 
reasons: in the bridge geometry they interact more with 
CO and this perturbs a lot the |m| = 2 PDOS on the 
neighbouring Pt; in the substitutional there is no CO or- 
bital with the matching symmetry and thus the |m| = 2 
channel is totally blocked by the molecule. Anyway, since 
these states fall below the Fermi level, their interaction 
with the molecule does not influence the conductance. 
The FR conductance does not differ substantially from 
the SR value because the \mj\ = 5/2 channels (which 
could in principle give rise to an increased conductance, 
since the corresponding FR band approaches the Fermi 
level) are almost completely blocked both in the upright 
bridge and in the substitutional case, and because there 
are no SO-induced gaps near the Fermi level. 

VI. CONCLUSIONS 

We have explored the bonding mechanism of a CO 
molecule onto a Pt monatomic nanowire. We found 
that in the bridge (and also in the on-top — i ^) config- 



uration the Blyholder model is appropriate to describe 
the electronic structure and the molecule-wire bonding. 
The HOMO and LUMO molecular orbitals are mainly 
involved in the bonding and strongly hybridize with Pt 
states forming bonding/antibonding states. As a con- 
sequence, from the 5(J molecular orbital some charge 
transfers to the wire, which in turn back-donates some 
of the charge, partially filling the 27r* orbitals of the 
molecule. The substitutional geometry has quite a high 
energy hence it could be a realistic configuration only 
when the wire is under high strain and almost broken. 
Also in this case there occurs a hybridization between 
the Pt states and the CO molecular orbitals which gives 
rise to a donation/back-donation process. We showed 
that the inclusion of SO coupling, which strongly mod- 
ifies the electronic band structure of the Pt wire, does 
not change very much the interaction mechanism and 
strength between the CO molecule and the wire. 

The adsorption of CO on bridge (and also on-top) ge- 
ometry should not affect much the chain ballistic con- 
ductance, since the Pt states mainly involved in the in- 
teraction with CO — and therefore partially or totally 
refiected — are located mainly below or above the Fermi 
level. In the substitutional geometry instead the molecule 
can partially transmit only |m| = 1 states, while all other 
channels are blocked. 

Our findings confirm the speculation by Kiguchi et al.— 
based on the measured conductances of transition metal 
nanocontacts, namely that the Blyholder model rules 
the adsorption strength of CO. Regarding the transport 
properties, our results agree with those of Strange et al.— 
in predicting that a considerable reduction of the ballis- 
tic conductance can be obtained when CO goes substi- 
tutional, and that the main conducting channel comes 
from a Pt \m\ = l-27r* hybridization (while m = states 
should not be conducting). Unlike previous calculations, 
we carefully examined SO effects on the intrinsic (tip- 
less) conductance and find that, for the specific case of 
the non-magnetic Pt nanocontact with an adsorbed CO 
molecule, a treatment of the electronic structure at the 
SR level was basically adequate to catch the effects on 
the conductance when bonding CO to the nanowire. 
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We used as reference SR (FR) all-electron configurations 
2s^2p^ (2s^/22p?/2) and 2s^2p'* (2s^/22p?/22p3/2) for C and 
O, respectively. The core radii for all channels are (1.3, 1.6) 
and (1.4, 1.6) in the C and O PP, respectively. The 3d chan- 
nel potential was treated as local with a core radius of (1.3) 
(C) and (1.4) (O), respectively. When two core radii are 
specified, that potential was pseudized in the US scheme, 
with the first radius representing the norm-conserving core 
radius and the second the US one. 
*® The all-electron configurations are the same used in the 
SR-LDA case. The core radu of Pt are (1.8, 2.2) for the 5d 
orbitals, (2.6) for the 6p and (2.4) for the 6s. For the 2s 
orbitals of C the core radius is (1.4, 1.6), while for all the 
other channels we used the same core radii of the LDA PP. 
The non-linear core correction^ has been included in all 
the PP. 

"^"^ We performed a test calculation with 25 Pt atoms for both 
bridge and substitutional SR geometries (see Sec. IHip and 
obtained very good agreement with the 17 Pt calculations. 
Around the Fermi energy the difference in the calculated 
transmission remains below 1%, while in the overall energy 
range the difference is at most 3%. 

In the structural optimizations we require that each com- 
ponent of the force on the C and O atoms gets lower than 
1 mRy/a.u. . In the bridge and on-top geometries the C 
and O atoms are constrained to move on the axis perpen- 
dicular to the wire [x axis) , while in the substitutional case 
we keep them aligned with the wire (hence on the z axis). 
In the tilted bridge configuration, the C and O atoms can 
move on the xz plane. 

In the SUBS-oxygen case we start with the O atom next to 
a Pt atom and the C atom distant from the Pt atom at 
the other side. For dpt-pt < 9.0 A the final configuration is 
the same as in the SUBS-carbon case (where we start with 
O distant from Pt), since the energy of that configuration 
is significantly lower. 

The atomic orbitals suited for the projection of states with 
symmetry can be chosen among spinors which are eigen- 
states of the total angular momentum , (with eigenvalue 
j{j -\- 1)) and of its projection along y, Jy (with eigenvalue 
jy). These states, labeled with {j,jy), transform according 
to the or irreducible representations. 



